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1. Introduction 

1.1. Background 

The Coulomb force is the dominating interaction in atomic and molecular physics. It 
is therefore the underlying force in chemistry and biology. Describing reactions with 
charged particles is thus an essential task for theoretical atomic and molecular physics 
as well as in chemistry. Despite this fundamental importance, solving the Schrodinger 
equation for even a three-body problem is a very difficult task. The asymptotic form of 
the wave function for three-body breakup is known to have very complicated structure. 
A recent review [Tj gave an overview of the most important articles which deal with 
the scattering problem of charged particles. The majority of these methods focus on 
solving the three-body problem without explicit knowledge of the three-body breakup 
asymptotics. Inspired by the methods of [2] and [21 E], in two recent studies we presented 
a new and rigorous method which can be used to solve the Coulomb scattering problem 
by using exterior complex scaling [5j E]. The former analyzes the single channel two- 
body problem while the latter indicates how the same formalism can be generalized 
to the full multi-channel three-body problem. In this contribution we demonstrate the 
construction of the formal part of the three-dimensional formulation of the three-body 
problem outlined in [6] by studying the three-dimensional two-body problem. 

1.2. Definition of the problem 

The scattering solution to the Schrodinger equation 

-A r + Vc(r)-A; 2 ]^ c (r,k) = (1) 
for the Coulomb potential Vc(r) = 2r]k/r has the well known form [TJ [SJ E] 

if>c(r, k) = r(l + ir/)e- 7r,?/2 e ir - k 1 F 1 (-ir/, 1, i(rife - r • k)). (2) 

Here the vectors r and k describe the position and the incident momentum. Their 
magnitudes are denoted by r and k. F and \F± are the Gamma function and the 
Confluent Hypergeometric function, respectively. Another way of representing the 
solution r/> c (r, k) is the partial wave decomposition [5] 

oo 

Vc(r,k) = — £(2£+ 1)^^(77, fcr)P,(cosfl). (3) 

fcr £=0 

Here cosd = r ■ k/rfc, on = argT(£ + 1 + vrj) represents the Coulomb phase shift, and 
Pi is the £th Legendre polynomial. The regular Coulomb wave function Fi(rj, kr) [10] 
obeys the partial wave equation 

/ d 2 £(£+1) 2r]k 



y (lr 2+^^ + ir- k2 ) F ^kr) = (4) 

and the boundary condition Fe(rj, 0) = at the origin. With this choice of coordinates, 
the wave function depends on the triad r,u = cos 9, k, such that ipc — ^Actj, u, k). 
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Although the explicit representations (12115]) have been known since the very earliest 
stages of Quantum Mechanics (7J |8] , the various approximations, which are based on the 
procedures necessary for suppression of the long range tail of the Coulomb potential, 
have been studied for many decades. The total number of publications on this subject 
is enormous, and so here we quote only those few which focus on the principal aspects 
of the problem [71 El [TTJ [121 EE] • Such procedures, called screening, are of substantial 
interest in view of their application to the scattering problem for more than two particles 
since in that case the exact solution of the Coulomb problem is not available [HI EE] . 

In principal, two kinds of screening procedures exist, i.e., the sharp and the soft. 
The sharp, which cuts off the Coulomb potential beyond some radius R, leads to the 
finite range potential 

V R (r) = V c (r)Q(R-r). (5) 

Here 9 is the standard Heaviside- function such that G(£) = 1(0), t > 0(t < 0). Soft 
screening methods imply multiplication of Vq by a smooth screening factor as, for 
instance, is done in the Yukawa potential W p (r) = Vc(r)e _r//p . The Coulomb wave 
function can be obtained in the limit p — > oo from the solutions of the Schrodinger 
equation with the potential W p . A regularizing factor is needed to obtain the correct 
limit. For the three dimensional Yukawa potential this factor is known analytically [13J. 

The situation with the potential Vr is more delicate. The representation for the 
Coulomb wave function through the solution for the screened potential Vr can easily be 
obtained at finite values of R for partial wave components [HI [111 [12] ■ The regularizing 
factor in this case is ^-dependent. The infinite sum over I should be computed to obtain 
the solution in three dimensions. An accurate analysis of the asymptotics of the partial 
wave series for the scattering amplitude as R — > oo was performed in [12] on the basis 
of distribution theory. That led the authors to the commonly accepted asymptotic 
regularizing factor e - l2r ? lo g 2fcR f or the scattering amplitude for the Vr potential. No 
extra term possessing singularities in the forward (backward) scattering direction were 
observed in [12], since the test functions used for the partial series summations were 
assumed to be vanishing in those directions. Actually, this requirement is not necessary. 
The problem of deriving the correct three dimensional expression for the wave function, 
which explicitly includes the Coulomb wave function ip c , was not the focus of the papers 
[HJ |12]. The paper [16] attempted to solve the problem in three dimensions directly 
by solving the three dimensional Lippmann-Schwinger equation with the potential Vr. 
However, the derivations made in [16] have been performed only for a particular value 
of the coordinate r = 0. It was not proven that the solution obtained in [16J is valid 
for all values of r. In the comment [TJJ it was shown by direct calculations that the 
three-dimensional result of [16] is erroneous. A further discussion of the results of [16] 
and their relation to the results of [12] can be found in |18j . 

The, to date, open situation of the cut-off Coulomb problem coupled with our own 
interests [5], related to the application of the complex rotation method for calculating 
the scattering states in the system of particles with long range interactions, stimulated 
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this research. The present paper is devoted to studying the scattering problem for a 
sharp cut-off Coulomb potential Vr and its complement 

V R = V c - V R . (6) 

In section two the partial wave equations are solved for a Vr potential by the 
conventional matching procedure at the point r = R and the infinite sum over angular 
momenta I is then evaluated in order to obtain the solution to the three dimensional 
Schrodinger equation. The main result of this paper is that we prove in the region r < R 
the wave function for the potential Vr in three dimensions is given by the action of an 
operator on the Coulomb wave function. This operator acts over the angular coordinate. 
The asymptotics of this operator is evaluated as R — > oo. Our derivation supports the 
form of the regularizing factor for the wave function which was proposed in [TT] for three 
dimensions but without a detailed proof. An extra term is found in the asymptotics of 
the scattering amplitude which possesses fast oscillations as a function of R and delta- 
functional singularity in the forward scattering direction. In section three we present the 
solution of the scattering problem for the Schrodinger equation with the V R potential. 
To the best of the authors' knowledge this represents the first time that this has been 
reported. This solution is used in section four to construct the three dimensional driven 
Schrodinger equation with the finite range potential Vr in the inhomogeneous term. As 
we demonstrated in our recent paper [5] , this equation is ideally suitable for applying the 
complex rotation method to solve the scattering problem with long range interactions. 



2. The scattering problem for the potential Vr 

The partial wave equation 

("^ + + Vn(r) - ^ v e (r, k) = (7) 

with the boundary condition ve(0, k) = and the asymptotics as r — > oo 

ve{r, k) ~ j e {kr) + A R£ hj(kr) (8) 

determines the scattering partial wave function for a given orbital momentum £. Here 
ji and h~l are the standard Riccati-Bessel and Riccati-Hankel functions |10j . The exact 
representation for ve(r, k) has a different form depending on whether the value of r is in 
or out the interval < r < R [TT]. For r G (0, R] one obtains 

v e (r, k) = ori F £ (r), kr). (9) 

For r > R the solution Ve takes the form 

v e (r, k) = j e (kr) + A Rl hj(kr). (10) 

At r = R both the function vi as well as its first derivative have to be continuous in r, 
i.e., d™ve(R — 0,k) = d^v^(R + 0, k), n — 0, 1. These conditions yield 

a R e = W R {j e , hj}/W R {F e , h+}, (11) 
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Am = W R {j £ , F e }/W R {Ft , h+}, (12) 

where W R {f,g} is the Wronskian f(r)g'(r) — f'(r)g(r) that is calculated at r = R. 
The phase shift 5m is then determined by the standard representation of the scattering 
amplitude A m 

A25 m _ "L 

A m = — j. — • (13) 
It is seen from (1111) and (1121) that the phase shift Sm can also be given by the argument 
of the amplitude am 

8 m = axga m . (14) 

Using the asymptotics of Riccati-Bessel functions as kR 3> £(£ + !) and the asymptotics 
of the regular Coulomb function as kR ^> £{£ + 1) + rf, one obtains the asymptotics of 
&m 

a m - ^ {ai ~ 1)log2kR) ■ (15) 

Therefore, the asymptotics of the phase shift 5m when kR 3> £(£ + 1) + rf reads [11] 

5m ~ o-e-rj log 2kR. (16) 

The above procedure describes how the partial waves vg(r, k) can then be constructed. 
Then the wave function v(r, u, k) is given by the infinite sum over momenta £ 

1 oo 

v{r, u,k) = — ]T(2£ + l)i%(r, k)P e (u). (17) 

fcr £=0 

This function satisfies the three dimensional Schrodinger equation ([T|) with the potential 
Vr taken instead of V c . 

Before proceed further we would like to point out that the convergence of the partial 
wave series for the scattering solutions should be considered with care, especially for the 
case of long range potentials [12J. The most reliable method is by using distribution 
theory. Consider an infinitely differentiable test function f(u) G C°°(— 1,1). By 
multiplying both sides of ( ITT)) with f(u) and integrating over u we obtain 



1 1 oo „i 

duv(r,u,k)f(u) = — ^2(2£ + l)i e v e (r,k) / duP,(u)f(u). 
-i krf^ J-i 

Introducing the Fourier coefficients with respect to the Legendre polynomials 



f £ =(2£+l) j\uP e (u)f(u), (19) 



equation ( ITS)) can be rewritten in the following form 



/i ^ 
duv(r,u,k)f(u) = —J2Mr,k)f e . (20) 
1 KT t=Q 
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The series on the right hand side is absolutely and uniformly convergent in the interval 
< r < oo, since the set of Fourier coefficients fe forms itself the absolutely convergent 
series and the terms of the series ( 1201) can be estimated as 

\vt(r,k)fe\<C\f t \, (21) 

where C is some constant. Thus, the leading term of the asymptotics of the series 
when r — > oo is now determined by the asymptotics of a certain number of coefficients 
Vi(r,k), whilst the tail of the series is negligible. The detailed description of using 
such an approach for the partial wave series summation can be found in [12]. In 
the following discussion we treat the partial series in the sense described above while 
assuming implicitly formulae such as (Tl8~l l20l) . We also extend this technique to the 
operators acting on the square integrable functions of the angular variable u. 



2.1. Properties of the solution for r < R 
For r < R the equation ( jTTl) yields 

oo 

v{r, u,k) = — + l^amFeiv, kr)P e {u). (22) 

Kr 1=0 

The right hand side of ( l22i) is the series in Legendre polynomials [in]. The polynomials 
Pi form an orthogonal and complete set of functions on the interval (—1,1) with respect 
to the scalar product 

(f\g} = J\ur(u)g(u), (23) 

where the asterisk indicates the complex conjugate. The orthogonality and completeness 
conditions for Pe are 

1 duPt(u)P x (u) = ^r^—Stx, (24) 



-1 



2£ + 1 



00 2C + 1 

£ —f-Pi(v)PM = 5 ( u ~ «')■ (25) 

£=0 z 

This set provides a basis for the C = L 2 (— 1, 1) space of square integrable functions 
on the interval (—1, 1), with (123]) as the inner product and with ||/|| = (flf) 1 ^ 2 as the 
norm. In the following derivation the elements of £ will be denoted as vectors, e.g. |/). 
With this notation ( 124"|) and ( 125|) take the abbreviated form 

{Pi\Px) = ^7^a, (26) 

E^HWd = I- (27) 
Here I denotes the unit operator in C. The equation (1221) now reads 

1 oo 

N r ^)> = 7rE(^ + i)iWifa»l^>- (28) 
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Here \v(r,k)) G C represents v(r,u,k) as a function of u. From the form of (j28p . this 
can be recast into 



\v(r,k)}= J2a m e-^ 2 -^\P e )(P e \x 

e=0 ' (29) 

i E(2A + l)i A e i ^F A (r/,fcr)|P A ). 

A=0 

By comparing this equation to (J3j) one identifies the right hand side of ( 1291) with the 
action of the operator 

00 2/ 4- 1 

aij = £a^e-^— — |P<)(P<| (30) 
£=0 2 

on the Coulomb wave function, which in £ is represented by the vector 

1 00 

\Mr,k)) = -^(2A + l)i A e^F A (r ? ,A;r)|P A ). (31) 

ftr A=0 

Thus, we have obtained the central focus of this part of the derivation, which establishes 
the relation between the solution of the Schrodinger equation with the sharply cut-off 
potential Vr and the Coulomb wave function for r < R. It has the form 

|tf(r,fc)) = a fl |Vc(r,A;)). (32) 

The inverse identity also holds true, yielding 

\Mr,k)} = aR 1 \v(r 1 k)}. (33) 

The exact form of a^, 1 is easily calculated from the orthogonality and completeness of 
Legendere polynomials 

2t+l, 



a 



R 



E fl «^-o-l p <)^l- ( 34 ) 



£=0 



The representation (iTTj) guarantees \or/\ 7^ and consequently the inverse is correctly 
defined by (1341) . The formulae (1321 1331) are the main results of this section, and they 
provide us with the representations of the Coulomb wave function in terms of the wave 
function v(r, u, k). These representations are valid for an arbitrary value of the screening 
radius R. 

A subsequent reduction of the complexity of the representations (J32l and ( J33l) can 
be observed if R — > 00 by studying the asymptotics of the operator a#. The operator 
a# can be simplified as R — > 00 with the help of the asymptotics (fl5j) of the coefficients 
am. The natural topology for calculating the asymptotics of a# is the strong operator 
topology in £, i.e. when the asymptotics of vectors a#|/) in C is considered as R — > 00. 
Construct a vector \gn) = a#|/) with an arbitrary |/) G C and then represent this 
vector as a sum of two terms \gn) = \giRi + \9r), where 

\gLR) = Ea m e- i °* 2 -^\P £ )(P e \f), 

^° (35) 
\9 L R }= E a^e-^|P)(P|/). 



£=L+1 
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The second sum can be made arbitrarily small by the choice of L. Indeed, due to the 
known properties of the Riccati-Hankel function and of the regular Coulomb function 
[10J it is seen that the coefficients ( fTTl) are bounded \am\ < C , then 



||^ f < C 2 g ^±^\(P e \f)\\ (36) 
e=L+i z 

The Parseval identity for the series in Legendre polynomials 

f*<2L^ mf) \* =W} (37) 

£=Q A 

shows that the series on the right hand side of (I3~T|) is convergent and, as a consequence, 
its tail 

i=L+l Z 

can be made arbitrarily small if L is taken to be sufficiently large. From this we conclude 
that the infinite sum on the right hand side of the inequality ( 1361) can be made arbitrarily 
small if L is large enough. Thus, for any small e > 0, there exists an integer L > 
such that for all L > L the inequality holds 

\\9rW 2 < C 2 e. (39) 

Consider now the vector \qlor)- Since L is finite there always exists a value of R such 
that the condition kR 3> L (Lq + 1) + rf is fulfilled. More precisely, the asymptotics 
of the Riccati-Hankel function and of the regular Coulomb function can be used for 
evaluating Wronskians in fTTTl) . as indicated in ffl5|) . to arrive at the inequality 

I a RLo - e i(CTL « -v i°g 2fcH) | < e i/2_ ( 40 ) 

Similar inequalities where L is replaced by I for all i < L are obviously also true. 
With these inequalities we get 

llbLo^-e-^^^J^II^ell/H 2 , (41) 

where 

L ° 2/ 4- 1 

Il =£^tP|P^|. (42) 

1=0 / 

Combining the inequalities f )39|) and f T4T|) together with the definition of \gn) we obtain 
the final estimate 

||(a* - e'^ 2kR I Lo )\f)\\ 2 < e(C 2 + ||/|| 2 ). (43) 

With this estimate, L can be extended up to infinity and the inequality with I instead 
of Il is also valid. Thus, the final result for the asymptotics of the operator is 
formulated as follows: 

Theorem 1. For any arbitrarily small e > there exists R such that the representation 
aR = e - ivlog2kR I + 0(e), (44) 
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holds true. Here the norm of a residual operator O(e) acting on any vector in C is of 
the order e as in . 

On applying Theorem 1 to the vector |/) = \ip c (r,k)) the asymptotic form of the 
relation fl32l) is obtained 



\v{r, k)) = e-^ 2kR \Mr, k)) + 0(e)|Vc(r, k)). (45) 

The alternative is given by 

\Mr, k)) = e^ 2kR [\v(r, k)) - 0(e)|Vc(r, fc)>] . (46) 

These two last formulae provide the strict basis for the problem of the asymptotic 
regularization [Til [T6] . which is needed when constructing the three dimensional 
Coulomb wave function from the solution of the Schrodinger equation with the screened 
Coulomb potential V R . 

Although the asymptotic regularizing factor e 1??log2A;R is numerical and does not 
depend on the angular variable u, the general factor sl r is the operator in the angular 
space C This operator connects the two solutions v(r, u, k) and V ; c( r , u -, k) to the same 
equations for r < R. This does not lead to a contradiction since the operator a#, or 
more precisely its extension on the three dimensional space for which we will keep the 
same notation elr, commutes with the Hamiltonian He = — A r + Vc(r) 

a R H c = H C 3L R . (47) 

Thus, if the function v(r,u, k) obeys the Schrodinger equation for r < R 

(H c - k 2 )v(r,u,k) = 0, (48) 

then one obtains 

(H c — k 2 )v(r, u, k) = (H c - k 2 )a R ip c (r,u, k) = 

= a R (H c -k 2 )Mr,u,k)=0. 1 } 

For the sake of completeness it is worth giving the explicit representation of the 
operator a^ acting on the functions of the angular variable 

& R f(r, u,k) = J ^ du a R (u, u')f(r, u, k). (50) 

Here the kernel a R (u,u') in accordance to (!30l) is represented by 

00 21 + 1 

a R {u,u') = J2a R ee~^^—P e (u)P e (u'). (51) 

This kernel is identical to the function (5) from [17] after respective unification of 
notations. 
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2.2. Properties of the solution for r > R 

For r > R the function vg(r, k) has the form ( JTU1) . The summation over i leads to the 
three dimensional solution 

v(r,u,k) = e irku + v sc (r,u,k) (52) 

where the scattered part has the form 

1 oo 

v sc {r,u,k) = -r- E( 2£ + l ) {tA ™ hj(kr)Pt(u). (53) 

fcr £=0 

Using standard arguments one arrives at the asymptotics of vsc as kr — » oo 

w sc (r, u, k) ~ fc)e ifcr /r (54) 

with the partial wave representation for the scattering amplitude Ar(u, k) given by 

1 oo 

A R {u, k) = - J2{2£ + l)A m P e (u). (55) 

If R — > oo the amplitude Ar(-u, fc) can be represented through the Coulomb scattering 
amplitude. In order to demonstrate this, one needs to use the representation of 
the amplitude through the phase shift ( ITBT) and its asymptotics ( |T6i) . The following 
representation is first derived 

1 00 pi2(o-^-77log2fc_R) _ 1 

M«> k )~Z E( 2 ^ + 1) 7p PtM- (56) 

By using the standard partial wave expansion of the Coulomb scattering amplitude A c 
Mu, k) = - £(2£ + lf——P e (u) (57) 

the representation ( 1561) can be transformed to the form 

ft) ~ e- i2r?los2/c/? A c (u, fc)- 
-|e- i "^ 2 ^sm(77log2A;i2) g (58) 

The sum in the second term of (158]) can easily be evaluated with the help of the 
completeness of the Legendre polynomials (}2"5]) and by taking into account that Pe(l) = 1 
00 2/ 4- 1 

E-^ L -W = *(!-«)■ ( 59 ) 

Here the delta-function is understood as in [12] 
■I 



du6(l -u)f{u) = /(l). (60) 

/-1 

Introducing (1591) into the formula (158|) we obtain the final form of the asymptotics of 
the amplitude (1551) 

fc) ~ e- 2i,?log2fcR Ac(M, fc) - l e - [rilog2kR sinolog 2^) - 1). (61) 

This is one of the main results of this section. One can recognize in the first term of 
(IB"Tj) the regularization factor, which was derived in [12] • The second term with the 
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strong delta function singularity was not known until now and, as demonstrated in the 
analysis shown above, its exact form can be obtained by the accurate summation of the 
all partial terms. 

3. The scattering problem for the potential V R 

In this section, the approach detailed in paper [5] is followed to construct the solution 
to the partial wave equation 

/ <l ~ iU j) +V R (r)-k 2 )w e (r,k) = (62) 



for the potential V R . The exact representation for We(r, k) is of the form 

wt(r, k) = afji(kr) (63) 

provided r < R. For r > R this gives 

wi(r, k) = e iae F e (r], kr) + Afu + (r], kr). (64) 

Here u + (r], kr) = e~ lcri (Gi + iFi) and Gi is the irregular Coulomb function [IU] . Similar 
to (JTT 1 fT2l) . the parameters af and Af should be found by matching the representations 
( 163]) and fl64|) for the function wg(kr) and for its first derivative at the point r = R. 
This leads to the expressions 

af = e^W R {F e , ut)/W R (j t , uj), (65) 

Af = ^Wr^^/WrCjzM)- (66) 
The phase shift 5f is introduced by the standard representation for the amplitude Af 

Af = e^-^. (67) 

For large values of R such that kR ^> £(£+!) +r] 2 the asymptotics of the regular Coulomb 
function and the Riccati-Bessel function can be used to get the following asymptotic 
representations for amplitudes (1651) and fl66l) 

af ~ e ir?log2fciJ , (68) 

i2{rjlo g 2kR-a l ) _ i 

Af~e i2 °'- -. (69) 

2i 

The solution to the three dimensional Scrodinger equation is given by the sum over 
the momenta i as 

1 oo 

w{r, u,k) = — + (r, k)P t (u). (70) 

kr 1=0 

As in the previous section this function takes the special forms on the intervals < r < R 
and r > R. The £ vectors will be used as above for formulating results. For < r < R 
the vector \w(r, k)) can be represented in terms of the vector |^'o( r ; k)), which represents 
the plane wave ipo( r i u i h) — e vrku , as follows 

\w(r,k)) = a R \Mr,k)). (71) 
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The operator a R is represented by 

a* = 5>?^p|J*><P,|. (72) 

The asymptotics of this operator as R — > 00 can be evaluated in the same way as in the 
previous section. The final result should again be understood in the sense of the strong 
operator topology in C. It reads 

This formula shows that, asymptotically, as R — > 00 

a*^ 1 . (74) 
For r > R the function w(r, u, k) is given by the expression 

w(r, u, k) = i[> c (r, u, k) + w S c{r, u, k), (75) 

where 

00 

w sc (r,u, k) = Y / (^+^ e Afut(7 1 ,kr)P e (u). (76) 
1=0 

Using standard arguments we obtain the asymptotics of w S c(r,u, k) as kr — > 00 

w sc (r, u, k) ~ A R (u, k)e l{kr ~^ 2kr) /r. (77) 

Here the partial wave decomposition for the amplitude is expressed as 

1 00 

A R (u, k) = - J2(2£ + l)A R P,{u). (78) 
fc 1=0 

If R — y 00 the asymptotics of the amplitude A R (u, k) should be combined with the 
Coulomb amplitude A c (u,k) in order to form the total amplitude of the outgoing 
spherical wave for the function ( |75l) . In this case we obtain 

A c (u, k) + A R (u, k)~y e irilos2kR sinolog 2*^) 5{1 - u). (79) 

k 

Comparing to the representation flBTj) for the amplitude Ar(u, k) one finds the relation 

A R ~ _ e 2iv\o g 2kR A ^ 

which holds for large values of R. This relation can be proved independently by 
comparing the partial wave series for the amplitudes A R (u, k) and Ar(u, k). 

The opposite limit as R — > has a certain interest for the case of the potential V R . 
It is obvious that V R — > Vq. The similar effect can be expected for the limit of the wave 
function w(r,u,k) — > ip c (r,u,k). The proof is based on the following asymptotics for 
the coefficients af and Af 

a? ~ e ^Q(2£+l)!!, (81) 
A R„^ t CtIL( kR) 2t+2 t (82) 
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Here Cf is the standard Coulomb normalization factor [10]. It is easy to see that 
asymptotically as R — > 

afj e {kr)~e ia *F e (r],kr) (83) 
for all < r < R. When R — ¥ the function u + (rj, kr) becomes singular as r — > R 

u + [r] , kr) oc r~ l . (84) 
At the same time the amplitude Af behaves as 

Af ocR 2 « +1 \ (85) 
Hence, for all R < r < oo one gets 

max \Afuj(r),kr)\ = 0(R e+2 ), (86) 

re[_R,oo) 

which shows that the term Afu^(t],kr) vanishes faster than the leading term 
e 1(7e Fg(rj , kr) = 0(r e+1 ) when r — >• R and R — > 0. Using these estimates in (164"]) one 
readily arrives at the statement 

w(r,u,k) ~ tp c (r,u,k) (87) 

when R — > 0. 

4. The driven Scrodinger equation for the Coulomb scattering problem 

In our recent study [5] we demonstrated that the inhomogeneous partial wave 
Schrodinger equation for the scattered part of the wave function with purely outgoing 
boundary conditions can be constructed and then successfully employed for solving 
the scattering problem for long range potentials. The key element of this approach is 
the solution of the partial wave Schrodinger equation for the potential V R . The three 
dimensional approach can now be formulated with the help of the solution from the 
preceding section. The Hamiltonian He is represented by 

H c = -A r + V R (r) + V R (r) (88) 
and the wave function ip c as 

ipeiji u, k) = w(r, u, k) + ipnir, u, k). (89) 
Here w(r, u, k) is the wave function for the V R potential 

(— A r + V R (r) - k 2 )w(r, u, k) = 0, (90) 

constructed in the preceding section. The function ipnir, u, k) obeys the inhomogeneous 
(driven) equation 

(-A r + Vc(r) - k 2 )Mr, u, k) = -V R (r)w(r, u, k) (91) 
and the purely outgoing boundary conditions as r — > oo, r > R 

i/) R (r, u, k) ~ A R (u, k)^ kr -^ 2kr ^ /r. (92) 
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It is seen from the definition that the amplitude A R is given in terms of A R through 
(EE} by 

A R = -A R . (93) 

For r < R it is useful to employ the interpretation of the functions as vectors in C. In 
this notations ip R takes the form 

\Mr,k)) = \Mr,k))-a R \Mr,k)). (94) 
and the equation ( !9~T1) becomes 

(H c - k 2 )\Mr,k)) = -a R V R \Mr,k)). (95) 
Multiplying (|94"|) by (a 72 )" 1 and using fl33l) . one arrives at the representation 

(a*)" 1 ^*:)) = ( aRa R )- l \v(r,k)) - \Mr,k)h 

which can be reduced using the asymptotic relation (1741) for large values of R to 

(a R )-Vi?M)> ~ \ v (r,k)) - \Mr,k)). (96) 

By its construction, the function 

\Mr,k)) = (a R y 1 \Mr,k)} } (97) 

obeys the equation 

(H c - k 2 )\Mr, k)) = -V R \Mr, k)). (98) 

Following (I9"61) . the function <p R (r, k) for r < R and large values of R can be represented 
as 

\<f> R {r,k))~\v(r,k))-\iP Q {r,k)). (99) 

Equation (I9"8l is the desired three-dimensional driven equation, which can be used 
for solving the Coulomb scattering problem by the complex rotation method. The 
formulation of the scattering problem on the basis of the equation (I98I) obeys two 
necessary conditions, which are needed for the application of the complex rotation 
method, i.e. i) the solution (p R (r,u, k) has the purely outgoing asymptotics 

(f) R (r, It, k) ~ (sl^Ar e Wr-vlog2kr) / r ^qq) 

and ii) the inhomogeneous term in the right hand side of (1981) vanishes outside of the 
sphere of the radius R. 



5. Conclusion 

New results on the structure of the solutions to the three dimensional Schrodinger 
equation for the sharply cut-off Coulomb potential have been derived. For the potential 
V R , which coincides with the Coulomb potential for all r < R, it was found that the wave 
function is proportional to the Coulomb wave function up to an operator factor. This 
operator acts as an integral operator over the spherical angular variable. The operator 
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a# is reduced to the multiplication by the constant e - l2r i lo & 2kR only asymptotically as 
R — > oo. This result clarifies the domain of validity for unjustified assumptions about the 
proportionality factor, which was taken as a constant in [UJ and |16j . The asymptotic 
representation of the scattering amplitude for the Vr potential in the case where R — > oo 
in addition to the standard term e~ l2r,log2kR Ac also contains the extra term ( 16"TT) . It has 
fast oscillations as R — > oo and a strong delta-functional singularity in the forward 
scattering direction. To the best of our knowledge this formula has been derived here 
for the first time. The representation for the scattering amplitude obtained in the 
recent paper [16] was derived from the incorrect form of the wave function in the region 
r < R [T7j and cannot be considered as a contra result. The complete solution for the 
scattering problem for the potential V R is given in this paper for the first time. The 
formula ( l80l) supports the complementary character of the two potentials in the sense 
that V R — Vq — Vr and should be considered as the fact of the self consistency of our 
treatment. The three dimensional formulation of the driven Schrodinger equation, given 
in section four, opens the way for forthcoming applications in the three body systems 
along the line given in [6]. 
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